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Résumé — On propose une commande de systèmes multi- 
variables, de dimension finie, linéaires ou non, sans en 
connaître le modèle mathématique. Nos deux outils essen- 
tiels sont l'algèbre différentielle, et une estimation des déri- 
vées de signaux bruités, récemment mise au point. Les simu- 
lations numériques de deux exemples, l'un linéaire, l'autre 
non, valident notre démarche. 

Mots-clés — Systèmes linéaires multivariables, systèmes non 
linéaires multivariables, identification boîte noire, estima- 
tion, dérivées de signaux bruités, algèbre différentielle, cal- 
cul opérationnel. 

I. Introduction 

Cette communication présente une généralisation multi- 
variable de [I^, où étaient jetées les bases d'une commande 
sans modèle pour systèmes monovariables de dimension fi- 
nie, linéaires ou non. Est-il besoin de rappeler {cf. [29j) 
la difficulté redoutable d'obtenir un modèle mathématique 
fiable dans bien des situations quelque peu complexes ? 
Ainsi s'explique l'incroyable popularité industrielle des cor- 
recteurs PID, en dépit de réglages souvent malaisés {cf. [2], 
[S]). Notre démarche se distingue par sa méthodologie et sa 
« philosophie » des identifications de type « boîte noire » , 
telles qu'on les trouve dans la littérature (voir, par exemple, 
[37] et f26]). Elle repose sur une nouvelle approché, qui a 
permis l'estimation en temps réel des dérivées de signaux 
bruités 



Nous substituons aux modèles mathématiques décrivant 
les machines dans une plage de fonctionnement aussi large 
que possible des équations différentielles « phénoménolo- 
giques », valides sur un court laps de temps, et actualisées 
pas à pas. Nous écrivons un système multivariable, de di- 
mension finie, à m entrées u = {ui, . . . ,Um) et p sorties 
y = (yi, . . . , î/p), sous la forme 



„("i) 



Fl + CH,1M1 +■■■-!- Ol.mMm + /3l 



^^"p) = _|_ ap^iui H h ap,m«m + /3p 



(1) 



^ Cette approche, née en [2Ï] à propos d'identific atio n paramétrique 
linéaire en boucle fermée, a été étendue au s i gna l [20| . Elle conduit à 
un c han gement radical de paradigme {cf. [ÏÔ], [H]). Voir, par exemple, 
0, |23| pour des applications concrètes. 

^Voir [Ï2], [13, [ni, [Ï5], [22], [24], ^ pour les déjà nombreuses 
applications en automatique non linéaire et traitement du signal. 



- rij >l,j = 1, 

- aj^i, (3j G M, î 



,p, et, le plus souvent, iij = 1, ou 2 : 



— 1, . . . ,TO, j = 1, . . . ,p, sont des pa- 
ramètres constants non physiques, choisis par le prati- 
cien ; 



les Fj sont déterminés grâce à la connaissance de y 



(ij, i = 1, . . . ,771, j = 1, . 



- Si p ^ 777, on garde seulement m sorties de manière à 
obtenir un système carré, inversible. 

Le comportement désiré s'obtient par correcteurs de type 
proportionnel intégral généralisé (GPI) [19J autour d'une 
trajectoire de référence. Dans les deux exemples ci-dessous, 
on utilise, comme en [Î6], un PID ou un PI. 

Le §[ni rappelle comment obtenir, grâce à l'algèbre diffé- 
rentielle, les équations différentielles entrée-sortie d'un sys- 
tème non Hnéaire. Le §|ïnisur l'estimation des dérivées d'un 
signal bruité inclut certains des progrès récents sur la mise 
en œuvre numérique. On évoque au § IIVI les principes es- 
sentiels de notre identification boîte noire. Le S IVl contient 
deux exemples, l'un linéaire, l'autre non, et leurs simula- 
tions numériques. Une brève conclusion évoque quelques 
perspectives futures. 

II. Rappels sur les systèmes non linéaires 

A. Corps différentiels 

Un corps différentiel R est un corps commutatif, muni 
d'une dérivation ^, c'est-à-dire une application R ^ A 
telle que, V a,b eÀ, 

- -§^{a + b) = à + b, 

- -§1(0^) = àb + ai). 

Une constante c G est un élément tel que c = 0. L'en- 
semble des constantes est le sous-corps des constantes. 

Une extension de corps différentiels consiste en la 
donnée de deux corps différentiels £, telles que : 

- Ç £, 

- la dérivation de ^ est la restriction à de celle de £. 
Notons ^{S), 5 C £, le sous-corps différentiel de £ engen- 
dré par K et S. Supposons £,/Â finiment engendré, c'est- 
à-dire £ = .^('S'}, où S est fini. Un élément ^ G £ est dit 

^Voir [5], |27j pour plus de détails et, en particulier, [5] pour des 
rappels sur les corps usuels, c'est-à-dire non différentiels. Tous les 
corps considérés ici sont de caractéristique nulle. 



différentiellement algébrique par rapport à Â si, et seule- 
ment si, il satisfait une équation différentielle algébrique 
• ■ • j C^"'') = 0, où P est un polynôme sur ^ en n + 1 
indéterminées. L'extension 2/Â est dite différentiellement 
algébrique si, et seulement si, tout élément de £ de différen- 
tiellement algébrique par rapport à Le résultat suivant 
est important : Z/^ est différentiellement algébrique si, et 
seulement si, son degré de transcendance est fini. 

Un élément de £ non différentiellement algébrique par 
rapport à ^ est dit différentiellement transcendant par rap- 
port à A. Une extension 2,/^ non différentiellement al- 
gébrique est dit différentiellement transcendante. Un en- 
semble {^t G £ I i e 3} est dit différentiellement algé- 
briquement indépendant par rapport à A si, et seulement 
si, aucune relation différentielle non triviale par rapport 
à Â n'existe : QC^f"'"*) = 0, où Q est un polynôme sur 

implique Q = 0. Deux ensembles maximaux d'éléments 
différentiellement algébriquement indépendants ont même 
cardinalité, c'est-à-dire même nombre d'éléments : c'est le 
degré de transcendance différentielle de l'extension £,/A. Un 
tel ensemble est une base de transcendance différentielle. 
Enfin, £,/A est différentiellement algébrique si, et seulement 
si, son degré de transcendance différentielle est nulle. 

B. Systèmes non linéaires 

Donnons-nous un corps différentiel de base k. Un sys- 
tèm^ est une extension différentiellement transcendante 
de K/k, finiment engendrée. Soit m son degré de trans- 
cendance différentielle. Un ensemble de commandes (in- 
dépendantes) u — est une base de trans- 
cendance différentielle de K/k. L'extension K/k{u) est 
donc différentiellement algébrique. Un ensemble de sorties 
y = (î/i, . . . , î/p) est un sous-ensemble de K. 

Soit X = {xi,...,Xn) une base de transcendance de 
K/k{u), de degré de transcendance n. Il en découle la re- 
présentation d'état généralisée : 

A,{x,,x,u, . . . .u'")) = 
Bf,(yK,x,u,...,uW) = 

OÙ L — 1, . . . , n, K — 1, • ■ ■ sont des polynômes 
sur k. 

La représentation entrée-sortie suivante résulte du fait 
que j/i , . . . , j/p sont différentiellement algébriques par rap- 
port à k{u) : 

$,(y,...,y(^.),u,...,u(*.)) = (2) 

OÙ ^j,j — l,...,p, est un polynôme sur k. 
Venons-en à l'inversibilité entrée-sortie : 

- Le système est dit inversible à gauche si, et seulement 
si, l'extension /c(u,y)/fc(y) est différentiellement al- 
gébrique. C'est dire que l'on peut récupérer l'entrée 
à partir de la sortie grâce à un système différentiel. 
Alors, m < p. 

- Il est dit inversible à droite si, et seulement si, le degré 
de transcendance différentielle de k{y)/k vaut p. C'est 
dire que les sorties sont différentiellement algébrique- 
ment indépendantes par rapport à k. Alors, p < m. 

*Pour plus de détails, voir [6], [Ï7], [35j . [36| . Rappelons que l'on 
emploie aussi l'algèbre difFérentielle dans diverses que stio ns relatives 
à l'identification et l'observabilité (voir, par exemple, [28]). 



Le système est dit carré si, et seulement si, m = p. Alors, 
inversibilités à gauche et à droite coïncident. On dit, si ces 
propriétés sont vérifiées, que le système est inversible. 

Remarque 1 Supposons notre système inversible à gauche. Le com- 
portement qualitatif de ^Ej), considéré comme système d'équations dif- 
férentielles en u, y étant donné, permet de définir le déphasage non 
minimal (voir aussi \25^). 

III. Estimation des dérivées temporelles 
À la série convergente x{t) = X]„>oa„^, a„ S 
C, correspond la série opérationnelle converg ent(ll X — 
Sn>o ï^TT- Avec le développement de Taylor tronqué 

XN{t) Yl,n=o °^nïïT' l'équation différentielle ^wtt ^ xn = Q 
correspond dans le domaine opérationnel à 

s^+^a'jv - s^xn{0) - s^-^iiv(O) ... - a:^'(0) = 

Les dérivées à l'origine x^^^ [Q) sont ainsi obtenues à partir 
du système d'équations linéaires 

{4"" (0) + 4"'"'' (0)- + ■ ■ ■ + -ivWs'^} = 

m — 0,...,N, V ^ + 1. Ce système étant triangu- 
laire avec des éléments diagonaux non nuls, les paramètres 
a;^''(0), et, par conséquent, les coefficients oo, . . . , cat sont 
linéairement identifiables [Î3], [2Ï]. Remplaçons xn par 
X dans ^ : on obtient ainsi V estimée opérationnelle 
[a:W(0)]e„ de :rW(0). 

Pour le passage au numérique, il suffit, selon les règles 
usuelles du calcul opérationnel {cf. [3Ô], [H], [33]), de rem- 
placer en ^ 

- pr, a > 1, c e C, par c^^^, t > ; 

- ^f-f par l'intégrale itérée d'ordre a 

/ / ■■■/ (-l)"T"a;(T)(it„_i---dtidT = 
Jo Jo Jo 

/ (t-r)-V"x(r)dr (4) 
(« - 1)! Jo 

Notons, [a;*^*^ (0)]e„ (i) V estimée numérique ainsi obtenue de 
a;'^*-'(0), pour un temps d'estimation t. La mise en œuvre 
repose sur le résultat suivant : 

lim[xW(0)]e„W = „lim [xW(0)],„(i) = xW(0) 



Remarque 2 Les itérations des intégrales produisent une moyennisa- 
tion, donc un filtrage passe-bas, qui permet d'atténuer les bruits (voir 



Remarque 3 La fenêtre temporelle d'estimation peut être choisie très 
petite, ce qui permet une implémentation en temps réel. 

^Voir |3Ô], [3Ï]. Renvoyons à [Ï3], [H] pour plus de détails. 



IV. Procédures d'identification boîte noire 

On requiert les propriétés suivantes pour éviter, notam- 
ment, toute boucle algébrique : 

1. On suppose le système inversible à gauche. Si le nombre 
de sorties est strictement supérieur à celui des entrées, 
c'est-à-dire p ^ m, on choisit m sorties pour obtenir un 
système carré inversible, et construire 

2. L'ordre de dérivation nj en |(T]) est relié à ^ par nj < Nj 



et, plus précisément, par 



dy'' 



^0. 



3. La valeur numérique de Fj, égale à y 



f3j est obtenue grâce à la discrétisation 



E 



ajAUi{K - 1) - Pj 



OÙ [•(K)]e désigne l'estimée au temps n. 
Le praticien suit les étapes suivantes : 

1. choix de m sorties si p > m ; 

2. choix des rij en (HJ ; 

3. choix de trajectoires de référence pour les yj, ainsi qu'il 
est usuel en commande par platitude (voir, par exemple, 

m, m, [310; 

4. choix des paramètres aj^i en ^ afin que que la grandeur 
des commandes u soient convenables ; 

5. choix des paramètres fij en ([T]) différent de si u n'ap- 
paraît pas Hnéairement en ([2]), c'est-à-dire |^(u = 0) = 0. 

Remarque 4 Avec des systèmes à déphasage non minimal, nos procé- 
dures peuvent conduire à des valeurs divergentes des Uj pour t grand, 
et, donc, à des valeurs numériquement inadmissibles des Fj . 

V. Deux exemples 

On utilise les modèles mathématiques ci-dessous pour les 
besoins évidents des simulations numériques. 

A. Système linéaire 

Soit le système linéaire à deux entrées et deux sorties, 
avec pôles instables et large spectre : 



î'I ~ (s+0.01){s+0.1)(s~l)3'"l 

~ (s+0.003){s-0.03){s+0.3)(s+3) "1 

(3+0.004)(s+0.04){3-0.4){3+4) "2 

Après quelques essais, nous choisissons ([T]) sous la forme 
découplée : 

yi = Fi + lOîti m = F2 + 10u2 

La stabilisation autour d'une traje ctoir e de référence est 
assurée par un régulateur PID (voir [16]) : 



(5) 



"1 = II) (yî - ^1 + Kpiei + Kii fei+ Koiéi) 
■"2 = ^ {§2 - F2 + Kp2e2 + K12 / 62 + Ko2é2) 

où 

- Kpi = 1, Kn = Kdi = 0, Kp2 = K12 = 50, Kd2 = 10 ; 

- î/*, 2/2 ^oiit les trajectoires de référence; 

- ei = î/î - yi, 62 = y2 - y2- 



^C 'est donc une commande prédictive sans modèle. Renvoyons à 
|18j et [7] pour les avantages de la platitude en prédictif avec modèle. 



Le comportement en suivi de trajectoires, avec bruit de sor- 
tie additif (loi normale Af(0,0.01)), est boijï|. Les estima- 
tions des signaux nécessaires à la synthèse de la commande 
sont présentées dans les figures [2]-(b) à[2]-(h). Notons la dif- 
férence d'échelle entre ui (figure [2]- (g)) et U2 (figure [2l- (h)), 
d'où l'inégalité des paramètres des correcteurs |[5]). 

La figure [3l-(b) montre le comportement du système 
en appliquant une commande PID plus « traditionnelle » , 
c'est-à-dire en posant Fi = F2 = 0. La comparaison des 
figures [3]- (a) et[3]-(b) est éloquente. 

B. Système non linéaire des trois cuves 

Le système des trois cuves de la figure [U très populaire 
dans la communauté du diagnostic!!, vérifie les équations : 



Xl 

±2 

±3 

y\ ■- 
y2 ■■ 
ys 



-Cisign(xi - X3)^\xi - xs\ + ui/S 
C3Sign(a;3 - X2)^/\x3 - X2\ 
-C2Sign(x2)^|x2 | -I- U2/S 
Cisign(a;i - Xi)y/\x[ 



2^3 1 



-C3sign(x3 - X2)y/\x3, - X2\ 



XI 
X2 
X3 



C„ = (l/S)./i„.SpV29,ra= 1,2,3 ; 

S = 0.0154 m (section des cuves) ; 

Sp = 5.10"^ m (section des tuyaux inter-cuves) ; 

g = 9.81 m.s~'^ (accélération de la pesanteur) ; 

Ml = ^^3 = 0-5, M2 = 0.675 (coefBcients de viscosité). 

Selon les recommandations du § \W\ on construit |(T]) , dé- 
couplé comme au S IV-AI : i/i = Fi-|-200wi, î = 1, 2. La figure 
|4l-(a) fournit le suivi de trajectoires. L'estimation des déri- 
vées (figure [4l-(b)) possède un comportement remarquable 
en dépit du bruit additif de mesure, de même caractéris- 
tique qu'au S IV- Al Les commandes nominales (figure |4l-(c)) 
sont assez proches de celles que nous aurions calculées en 
utilisant la platitude (voir |15j). Elles sont complétées par 
des correcteurs PI 



1 

200 



lOei +2.10" 



1,2 



OÙ y* est la trajectoire de référence, Ci = y* — yi. Pour 
évaluer nous utilisons yi débruité (voir figure|l]-(d)) selon 
les techniques du §[1111 (voir aussi [13j, ^14]). 

VI. Conclusion 

Les règles ardues d'identification du S IIVI seront précisées 
dans le futuiH. On exposera bientôt des résultats encoura- 
geants sur le déphasage non minimal ainsi que ceux sur 
V égalisation aveugle {cf. [9j), qui est, en un certain sens, le 
pendant en signal de la commande sans modèle. 

^ C'est pourquoi nous espérons, comme déjà dit en |16| . que nos 
méthodes pourraient fournir une alternative efficace à celles sur la 
réduction de modèles (voir, par exemple, [î]). 

^Renvoyons à [U] pour plus de détails et les ré féren ces bibliogra- 
phiques. Que l'on nous permette de rappeler que [H] propose, sans 
doute pour la première fois, le diagnostic, la commande et la reconfi- 
guration d'un système non linéaire à paramètres incertains. 

^EUes ne peuvent être toutes déduites de considérations purement 
mathématiques. L'expérience pratique y joue bien entendu, comme 
pour les PID, un rôle considérable. 




Fig. 1: Système non linéaire des trois cuves 



Une mathématisation élaborée, comme celle présen- 
tée ici, afin d'abandonner une modélisation aussi « glo- 
bale » que possible dans une discipline empirique, comme 
l'automatique, semble nouvellj^. Il conviendrait d'en ex- 
plorer les implications épistémologiques, ne serait-ce que 
pour les liens entre physique, complexité, commande et ré- 
solution temporelle {cf. |32j). 
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Fig. 4: Système non linéaire des trois cuves 
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